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Abstract 

We have analysed here the role of the geometric phase in dynamical mechanism of quantum 

phase transition in the transverse Ising model. We have investigated the system when it is driven 

at a fixed rate characterized by a quench time Tq across the critical point from a paramagnetic 

to ferromagnetic phase. Our argument is based on the fact that the spin fluctuation occurring 

during the critical slowing down causes random fluctuation in the ground state geometric phase 

at the critical regime. The correlation function of the random geometric phase determines the 

excitation probability of the quasiparticles, which are excited during the transition from the inital 

paramagnetic to the ferromagnetic phase. This helps us to evaluate the number density of the kinks 

_ 1 

formed during the transition, which is found to scale as ^ . In addition, we have also estimated 
the spin-spin correlation at criticality. 
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I. INTRODUCTION 



The celebrated geometric phase, commonly known as the Berry phase |l| arises when a 
quantum system is adiabatically transported around a closed circuit in a parameter space. In 
this transportation, besides the familiar dynamical phase, a geometric phase appears when 
due to nonintegrability some variables fail to return to their original values. The geometrical 
meaning of this phase is precisely the holonomy in a Hermitian line bundle and naturally 
defines a U{1) connection The Berry connection and the geometric phase can be defined 
even if the condition of adiabaticity is relaxed [3!]. Subsequently, the generalization of the 
Berry phase was introduced to the non-Abelian case where a collection of degenerate 
states were cycled. An analogous phase in classical mechanics was also pointed out js, 6|. 
The case with non- Hermitian evolution was also studied [t]. The interpretation ^| of the 
driving parameters as dynamical variables, whose evolution was influenced by the same 
geometrical objects as the Berry phase, has also been noted. This made clear |9| that the 
reaction of the geometric phase on the parameters takes the form of an abstract magnetic 



field which was explored earlier 



10|. A classic connection was made [11] between the Berry 



phase and the geometric phase in beams of light 



12| . which was formulated for two state 



systems in the optical picture by Pancharatnam ISj. The concept of the Berry phase has 
now become very important in quantum mechanics with wide applications in various fields 



14j . Il5| . In particular, geometric phases have been associated with a variety of condensed 



matter and solid state phenomena 16 



of experimental investigations also 



14 



22 L Besides, theoretical analyses it has been a subject 



23 



29|. 



A characteristic feature of all non-trivial geometrical evolutions is that it demonstrates 
the presence of non-analytic points in the energy spectrum where the state is degenerate. 
It may be mentioned here that before the discovery of the geometrical phase, the sign 



30, 



3lj. The 



change associated with the molecular electronic degeneracies was understood 
presence of degeneracy at some point is accompanied by the curvature in its immediate 
neighbourhood and a state evolving around a closed path can detect it. These aspects 
indicate an intuitive connection between the criticality around the degeneracy point and 



the geometric phase. Recently, the geometric phase has been studied |32|-|34| in the context 
of quantum phase transitions (QPT)[35|, |36| to investigate the regions of criticality. QPTs 



exhibit the crossover of ground states at a critical value of some external parameter at zero 



2 



temperature and are caused by the presence of degeneracy in the energy spectrum. The 
geometric phase, which measures the curvature associated with the degeneracy point, can 
contemplate the energy structure of the system and essentially captures certain characteristic 
features of QPT. Pachos and Carollo [37] have exhaustively studied the role of geometric 
phase as a tool to probe quantum phase transitions in many-body systems and suggested the 
feasibility of the experimental realization also. The geometric phase has also been used as 
a topological test to reveal quantum phase transition j34|. Recently it has been shown that 
the geometric phase associated with the ground state of the XY model exhibits universality 
in the critical properties and obeys the scaling behavior near the critical point [33], which 
established a significant connection between the geometric phase and the quantum phase 
transition. 

The inter-relation between the geometric phase and quantum phase transition has so far 
been concerned with the detection of criticality through this phase and focus on equilibrium 
scalings of various physical quantities near the critical point. However, till now, the role of 
the geometric phase in the dynamical mechanism of QPT has not been explored. It is our 
novel attempt here to show that the geometric phase plays a crucial role in the dynami- 
cal mechanism of the QPT when it is induced by a linear quench, where non-equilibrium 
dynamics while passing through a critical point becomes relevant. The Kibble-Zurek mech- 



anism (KZM) 38|, |39| which is a universal theory of the dynamics of the second order phase 



transition, has recently been explored in the zero temperature limit to study the dynamics of 
the quantum phase transition . The behaviour of the quantum Ising model in one dimension 
las been investigated during a quench caused by gradual tuning of the transverse bias field 
40 . 141 I. The system is driven at a fixed rate characterized by the quench time Tq across the 



critical point from paramagnetic to ferromagnetic phase. In agreement with KZM (which 
recognizes that evolution is approximately adiabatic far away, but becomes approximately 
impulse sufficiently near the critical point), quantum state of the system after the transition 
exhibits a characteristic correlation length ^ proportional to the square root of the quench 
time Tq-. ^ = ^Jjq. The inverse of this correlation length is known to determine the average 
density of defects formed during the transition. We focus here on the role of the geometric 
phase in the dynamical mechanism of QPT induced by a linear quench, and have established 
a nontrivial connection between KZ mechanism and the geometric phase. It is pointed out 
that the quench induced QPT in a spin chain generates spin fluctuation which causes the 
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randomization of the time dependent magnetic field so tliat the geometric phase acquires 
random values. The two point correlation of randomized geometric phase corresponds to 
the excitation probability of the quasiparticles. From this we have estimated the number 
density of defects generated in the transverse Ising model during the phase transition. 

The plan of the paper is as follows. In section 2, we review the geometric phase (Berry 
phase) of the ground state of the transverse Ising model as discussed in |32| . Section 3 deals 
with the discussion of the geometric phase during critical slowing down induced by a linear 
quench. In section 4 we estimate the number density of kinks from the two-point correlation 
of the Berry phase factor. Section 5 is devoted to discuss the spin-spin correlation at the 
critical point. 



II. GEOMETRIC PHASE OF THE GROUND STATE OF THE TRANSVERSE 
ISING MODEL 

The Hamiltonian of a one- dimensional transverse Ising model is given by 

+M 

H=-jC£<+^^t (1) 

-M 

with M = ^^Y^, where (odd) denotes the number of sites, A is the external field and ex's 
are the standard Pauli matrices. We assume periodic boundary condition. In the transverse 
Ising model, the geometric phase of the ground state is evaluated by applying a rotation 
of an angle a around the z-axis to each spin jssi, Is^]. A new set of Hamiltonians Ha is 
constructed from the Hamiltonian ([T]) as 

Ha = U{a) H U^{a) (2) 

where 

+M 

U{a) = Yl exp{ia(T]/2). (3) 

j=-M 

We use the standard Jordan- Wigner transformation, which maps the spins to one dimen- 
sional spinless fermions, 

^^={U<^^] (4) 



4 



and the Fourier transforms of the fermionic operator 



-27rj k 



with k = -M,... + M 



(5) 



such that the Hamiltonian Ha can be diagonahsed by transforming the fermionic operators 
in momentum space and then using the Bogohubov transformation. This yields 



where 



7 -A 2ia ■ 

Ck = dk CO?,— - id^_^e smy, 



is the fermionic annihilation operator in the momentum space, with 

2TTk ^ 



COS ^ 

COS 6*4. = 



and 



Aa 



A — COS 



2'nk 



+ sin 



2'Rk 
1^ 



(6) 



(7) 



(8) 



(9) 



the energy spectrum. 

The ground state \g > of the system is expressed as 

M 



>= ]^(cos ^|0 >k |0 >_fe -iexp(2ia) sin ^|1 >k |1 >-a 



\9 



(10) 



k=l 



where |0 >fc and |1 >fc are the vacuum and single fermionic excitation of the k-th momentum 
mode respectively. The geometric phase of the ground state, accumulated by varying the 
angle a from to vr, is found to be [32|, |37 | 

d 



T = -i 



9 



da 



g) da = 7r(l — cos6'fc) 



(11) 



fe>0 



III. GROUND-STATE GEOMETRIC PHASE IN A VARYING MAGNETIC 
FIELD AND CRITICALITY 

We start with the review of the geometric phase of a spin 1/2 system in a varying 
magnetic field. The Hamiltonian of a single spin 1/2 system in the presence of an external 
time dependent magnetic field is given by 

H{t) = ^-Bit).a (12) 
5 



where a = {ax,(Ty,az) are the Pauh matrices and B(t) = Bon{9(t)), n{9,t) being the 
unit vector with n = (sin 6* cos 0, sin ^ sin 0, cos ^). If the external magnetic field is varied 
adiabatically, the instantaneous spin states can be expressed in the cr^ basis as 



9 9 

tn;^> = cos-\ tz> + e"*' sm-\ 



9 9 
U;t> = sin-| t^> - e^-^cos-l 



(13) 



where | t^>) I iz> ^^^e the eigenstates of the cr^ operator. 

For a cyclic time evolution i.e. for B(T) = B(0), apart from the dynamical phase, the 
eigenstates acquire a geometric phase also and we can write, 

I tn(r)>=e^V^^|tn(o)> (14) 

where the dynamical phase S = BQ{t)dt and the geometric phase (GP) 7 = / A''' ■ dX; 
A is the set of control parameters and = i <^n | Va| tn> is the so called Berry connection. 
In our case, A = {9, 0) and a straight forward calculation shows that 

7t = -7; = ^(1 - cos 6') (15) 

It may be pointed out that though the eigenenergies depend on Bo(t), the eigenstates depend 
on n{9{t)) only. The GP is proportional to the solid angle subtended by B with respect to 
the degeneracy B = 0. It may also be noted that fiuctuation in the external magnetic field 
will obviously induce fiuctuation in the geometric phase (7^ or 74^) of the corresponding spin 



eigenstate through fiuctuation in cos 6* j42|. 



This formulation may be extended for a chain of spin | system with XY type of inter- 
actions, which exhibits QPT for a critical magnetic field. Subsequently, we can analyse the 
dynamics of the geometric phase when there is a gradual slowing down of the external mag- 
netic field. We consider that the quantum fiuctuation responsible for a QPT induces spin 



fiuctuation 43j which eventually generates random fiuctuation in the ground state Berry 
phase. In the critical regime, the Berry phase of the ground state of the relevant system 
undergoes a random fiuctuation which is assumed to be small for a slow transition and large 
enough for a fast transition. In the present communication, we have used this specific feature 
of the geometric phase of the ground state of the transverse Ising model during its critical 
slowing down and showed that the two-point correlation of the Berry phase factors may be 



used to determine the excitation probability of the quasiparticles and hence the density of 
defects formed during the transition. 

Let the system governed by the Hamiltonian (1) be initially in its ground state at the 
large values of A >> 1, a paramagnet with all spins polarized along the x-axis. Gradually, 
when A is ramped down to very small values A << 1, there may be two degenerate ground 
states representing ferromagnets either with all up spins or all down spins along the z-axis. 
For a large number of spins i.e. N ^ oo, the energy gap at the critical point A = Ac = 1 
tends to zero implying that the system is excited while passing (with a finite speed) through 
the critical point and the system settles down in the final state with A = 0. As a result the 
system in the final state represents a quantum superposition of states as 

l--unmut > 

with finite domains of spins pointing up or down. The domains are separated by kinks 
or antikinks, where the polarization of spin changes its orientation. During the excitation, 
spin flips occur at different domains of the chain so that in these regions there are two 
nearest neighbour spins with opposite orientations. The average density of kinks depends 
on the transition rate. For a slow transition the number of kinks present is small. For a 
fast transition, the orientation of individual spins become random, uncorrelated with their 
nearest neighbours and the density of kinks becomes large. We would like to explore the 
situation when the time dependent field X{t) driving the transition is given by 

X{t){t < 0) = -- (16) 

where the quench time Tg is an adjustable parameter. We may note that before the quench 
{-t > Tq), the system is in the symmetric (paramagnetic) phase. During the quench time Tq, 
it undergoes a transition towards the symmetry broken ferromagnetic phase and finally at 
the end of the quench (t = 0) it settles down. 

The total time dependent Hamiltonian of the system may be written as 

H{t)^Ho{t)+HfUp{t) (17) 

If jV pairs of quasiparticles arc excited with momentum k^., — /coi with (i = 1,2, ...J\f) at a 
certain quench time t — Tq, the time dependent spin fiip interaction is given by 

H,uAt) = -2 J E (18) 

feo=l 

7 



with 

(j""^ = I i><t I or I t><i I 

and ^ko{t) being a time dependent dimensionless parameter characterising the spin flip 
interaction strength. Here the summation includes the spins which get excited during the 
transition and finally are fiipped with opposite orientations. This term corresponds to the 



standard Landau-Zener Hamiltonian 



44j |. The first term in eqn. (17) is given by 



Hoit) = -JiYl E ^^^Ui + KtM) (19) 

n 4=1 

Nn being the number of spins between the fiipped spins in the nth domain and n is the 
number of domains in the whole system. This corresponds to the Hamiltonian of the sur- 
rounding Ising spin chain. We have neglected the coupling between the fiipped spins with 
the surrounding Ising chain. 

In the process of critical slowing down certain spins which are excited, undergo fiuctua- 
tions around the z- axis. In fact, the spin fiuctuation is instrumental in a metal- insulator 
transition which in zero temperature limit corresponds to QPT [43]. When a spin gets ex- 
cited, it passes through the process of the change of its alignment around the z-axis until at 
the end of the quench at t=0 it settles with opposite orientation. This process of the change 
of alignment of spins at the quench time can be depicted as a fiuctuation of the spin when 
magnetic moments appear with random orientations. 

It may be noted from eqn. (15) that during critical slowing down the fiuctuation of the 
angle 6 associated with spin fiuctuation, makes the Berry phase a fiuctuating one. Indeed, 
when a quasiparticle associated with a spin rotates around a closed path, the Berry phase, 
which is given by the holonomy corresponds to the number of magnetic fiux lines enclosed 
by the loop. As during the quench, the time dependent magnetic field is randomized due to 
the spin fiuctuation, the Berry phases acquired by the spin eigenstates take random values. 

To evaluate the Berry phase acquired by the spin eigenstate associated with the quasi- 
particle undergoing excitation near the critical point we note that at the quench time Tq, 
spins tend to align along the z-axis as the initial paramagnetic system (A >> 1) transits 
to a ferromagnetic one (A << 1) passing through the critical point (Ac =1)- In fact the 



critical Ising model corresponds to a free fermionic field theory 



45l | with spins aligned along 



the z-axis. To determine the Berry phase at the quench time ■jkoi'i'q) for a spin eigenstate 



8 



corresponding to a quasiparticle with momentum ko, we note that initially at A >> 1 the 
system is in the paramagnetic state with all the spins aligned along the x-axis. We consider 
that the time dependent magnetic field rotates around the x-axis with an angular velocity 
Uq. At time t = Tg {jq = the spin state with momentum mode k^, picks up a geometric 
phase 

7fcoK) = 7r(l - cos^fc„(rg)) (20) 

with Ok^ijq) = 7r/2 indicating that at t = Tg, the geometric phase ■jkoi'Tg) acquired by the 
spin eigenstate is vr. 

We consider that during critical slowing down A/" pairs of quasiparticles with momentum 
modes kQ., —k^- with (z = 1, 2, ...A/") are excited with probability pj^.^ . The corresponding spins 
are fiipped with opposite orientations and settle down at t = implying 9koit = 0) = vr, such 
that the random variable 7fco(t = 0) takes the value 27r with probability p^iy In essence, the 
Berry phase acquired by a spin eigenstate at t = and at t = are given by the random 
values 

7M=7r, 7fc„(t = 0) = 27r (21) 

with probability p^o ^ind with with probability 1 — pk^- We can now estimate the number 
density of kinks generated in the system due to a quench induced QPT. 



IV. ESTIMATION OF THE NUMBER DENSITY OF DEFECTS 

For our further analysis, we define a (dimensionless) random variable 

Mt) = ^ (22) 
/vr 

such that at if: = Tg and t = 0, the Berry phase factors (pkoit) attain the values 

Kit = rg) = l/2, 0fc„(t = O) = l (23) 
Let us assume that the Berry phase factor (pkoit) follows the simplest stochastic differential 



equation 



d^ko (t) = -(^ko^ko it)dt + dri{t) (24) 

where UkQ is the frequency associated with the energy of the quasiparticle near the critical 
point and ri{t) denotes the fiuctuation. Let us consider ri{t) as the Gaussian white noise with 



moments 



< dri{t) > = (25) 
< dri{t)dri{t') > = 5{t-t')dt' (26) 

Using eqn.f l24p . and performing the averages over the noise, the correlations of the time 



dependent Berry phase factors can be evaluated 



47|. This gives 



<0fe„(t) > = (27) 
<K{t)Kit')> = ^e--'=o(*-*') (28) 

This leads us to derive the two-point correlation function of the Berry phase factor 
during the critical slowing down from t = to t = as 

<0fco(^.)0fco(O)>=^e-"'=°^' (29) 

The correlation is large for fast transition (small Tq) and small for slow transition (large r^). 

As JV pairs of quasiparticles are excited during the critical slowing down with momentum 
modes kQ.,—kQ. {i = 1,2, ...A/") with probability pkg, we note from relation (23) that the 
random variable \2(f)ko{t = Tg)(f)ka{t = 0)| picks up the value 1 with probability p^o 
with probability 1 — ■ This suggests the excitation probability to be given by 

Pko = 2 (0fco(^g)</'feo(O)) 

Using eqn.(29), this can be explicitly written as 

p,, = 2 {(f>MKm = e-^-or. = (h = 1) (30) 

One should note that in the final state at t=0, the flipped spins correspond to the values 
of 9ko = vr, with cos9ko = —1- Now if we choose the lattice spacing a = ^ = 1, in the 
critical region, with A = 1, the constraint cos9ko = —1 suggests from eqn. (8) that 

Afco = 1 - cos ko (31) 

The expression for the energy then yields 

ko 

efcp = 2JAfc(, = 2J(1 — cos/co) = 4Jsin^ — ~ Jk^ for small kg (32) 



10 



From this, we find 



Pfco = e-'^ ^ 'o^^ ih=l) (33) 
The number of kinks can now be evaluated as 

2X=Y,Pko (34) 

ko 

For large number of spins, i.e. in the thermodynamic limit as N ^ oo, the expectation 
value of the density of kinks is given by 

2A/' 1 1 1 

n =N^oo ^ = ^ y Pkodko = ^(2 J T^r^ (35) 

_i 

which shows that the density of kinks scales as \ In this novel framework, the dynamics 
of the geometric phase associated with the flipped spins helped us to estimate the density 
of defects produced in a system due to QPT and the result is found to be in concordance 
with that derived by a standard technique 41 1. 



V. SPIN -SPIN CORRELATION AT THE CRITICAL POINT 

We can now derive the spin-spin correlation at the critical point in terms of the space 
dependent correlation of the geometric phase. In the initial ground state, with A — ?■ oo, due 
to symmetry < cr^ >= and the transverse magnetization < cr^ >= 1. For a slow quench, 
the final magnetization, < a^' >— ?■ 0, and in the final state at A = 0, < o"^' >= 0. It may be 
noted that for a certain quench time Tg, < > will depend on the density of kinks . In the 
region between a kink-antikink pair, the direction of the spin orientation remains the same 
corresponding to the direction of orientation of the flipped spin. If the kink-antikink pair 
maintains more or less same distance from each other, the system represents a kink-antikink 
chain with a lattice constant ~ where ^ = ^/r^ is the Kibble-Zurek correlation length. 



Indeed the expected value < erf > for a particular quench time is [48 1 



<^^>-V^7^ (J=^) (36) 

which is valid for Tg » 1. We have just analysed how the variation in the direction of the 
spin orientation during the quench time Tg is reflected in the geometric phase factor (pk^^Tg). 
In fact, the direction of the spin orientation at the spatial position r denoted by site i near 

11 



the critical point determines the geometric phase 7(r) acquired by the spin eigenstate after 
a cychc evolution. If r' denotes the spatial coordinate of the spin at site i + R, then we can 
define the two-point correlation of the Berry phase factor near the critical point as 

=< <jy{r) <j,{r') > (37) 

where 0(r) = ^7(^)- The correlation f l37p effectively corresponds to the the spin-spin cor- 
relation at the critical point. As already mentioned, the Berry phase undergoes a stochastic 
fluctuation near the critical point. In terms of spatial variable, we consider that the Berry 
phase factor corresponding to the spin state at position r follows the stochastic differential 
equation 

d(f){r) = -^^dr + dri{r) (38) 

Here //(r) is the fluctuation and ^ is the KZ correlation length which is the relevant length 
scale in the system. We consider that ri{r) represents Gaussian distribution with white noise 
and satisfies the moments 

< dr]{r) > = (39) 

< dri{r)dri{r') > = 6{r — r')dr' (40) 

I — I 

From this the correlation function of the Berry phase factor can be evaluated [47] and is 
given by 



< (f){r) > = (41) 
1 

-f 
2 



< 0(r)0(r') > = -e~^^ (42) 



This implies that 



=< 0(r)0(r') >=le f = 0.5e f {h = J = 1) (43) 



This can be compared with the numerical estimate of spin -spin correlation at the critical 
point in ref. [48] 

^ 0M^-2mf ^^^^ 

which is valid for Tg » 1 and i? >> ^ 

It is noted that in contrast with the result (43), the relation (42) is a generalised one and 
does not involve any restriction on the values of Tq and R. It is observed that for very small 
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Tq {Tq << 1) the Spin -spin correlation tends to zero which is consistent with the fact that 
for fast transition the spins behave almost uncorrelated. Moreover, contrary to eqn. (43) 
we have found that there is a finite correlation for large Tq and large R. 



VI. CONCLUSION 

The KZ mechanism in the second order phase transition suggests that the formation of 
defects is driven by thermal fluctuation. In the analysis of spin models, the KZ theory 
has been extended to study quench induced quantum phase transition driven by quantum 
fluctuation. Equivalently, we have argued that when the KZ mechanism is encompassed 
in the quantum phase transition of spin systems, the defect formation is driven by the 
fluctuation of the Berry phase factor associated with the ground state. Recently, we have 
studied a XY spin chain in a slowly varying time dependent magnetic field and investigated 
the behaviour of the geometric phase during a linear quench caused by a gradual turning 
off of the magnetic field 49|]. In the present formulation, we have argued that in quantum 
phase transition induced by a quench the Berry phase factor undergoes a random fluctuation 
at criticality. The two-point correlation function of the Berry phase factor determines the 

dynamics of the QPT induced by a quench and controls the formation of defects. It is 

_ 1 

found that the number of kinks generated in this transition scales as Tq ^ where Tq is the 
quench time. The result is in agreement with that obtained through other formulations. 
Besides, the space dependent correlation of the random Berry phase gives us the estimation 
of the spin-spin correlation at the critical point. It is hoped that this work may initiate 
investigations in relation to the role of the geometric phase in the dynamical mechanism of 
quantum phase transition. 
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